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Abstract 

We re-examine soft dilepton emission from a weakly coupled Quark-Gluon Plasma. We show 
that Braaten, Pisarski, and Yuan's result that the dilepton rate rises as 1/% (and the spectral 
weight scales as l/q°) at small energy q° <C gT is correct, but that the coefficient they found for 
this behavior is not correct, because their analysis was incomplete. At still smaller scales, the 
behavior moderates to oc for q° < g 4 T, consistent with a finite electrical conductivity. We 
evaluate the spectral weight in the q° ~ g 4 T region by kinetic theory techniques and show that it 
satisfies a sum rule, which makes the determination of electrical conductivity from the Euclidean 
correlation function very challenging. 
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I. INTRODUCTION 



It is believed that experiments at RHIC are producing the Quark-Gluon Plasma, which 
will be further probed by the heavy ion program at the LHC However, one of the com- 
plications in extracting the properties of the Quark-Gluon Plasma is that thermalization 
and rescattering tends to erase information about the state of the plasma at early times. 
Therefore, there has been an effort to understand so-called "hard probes," particles which 
are produced early during a heavy ion collision and either escape without further interaction, 
or interact with the high density plasma much more strongly than with the hadrons present 
at late times. Jets and heavy quarks are examples of the latter; photons and dileptons are 
examples of the former. For a comprehensive review, see for instance [1]. 

The thermal production of dileptons has long been viewed as a promising means to 
identify and study the Quark-Gluon Plasma [2, 3]. Theoretically, our two best tools for 
establishing the dilepton production rate by the plasma are perturbation theory and lattice 
gauge theory. Each is rigorous within some regime. Both suffer from limitations. In partic- 
ular, while perturbative tools are rigorous in the limit of small coupling (for a fixed energy), 
they probably do not show good convergence at the value of coupling constants relevant in 
realistic experiments. Furthermore, the application of perturbation theory has proven much 
more complicated than naively expected. In particular, when momenta are small (or close 
to the light cone [4]), the perturbative expansion requires a resummation, the Hard Thermal 
Loop (HTL) resummation, pioneered by Braaten and Pisarski [5, 6]. 

The thermal dilepton production rate appears to be under control for large, deeply time- 
like momenta. Perturbation theory works well; the leading and 0(a s ) subleading behavior 
is known, and there are no singularities encountered [7]. Further, lattice gauge theory tech- 
niques are applicable and show good agreement with perturbation theory [8]. However, the 
situation is less clear for small or nearly lightlike momenta. For large but nearly lightlike 
momenta, perturbation theory can be applied, but it requires rather extensive resummation 
[9]. It may be quite difficult to compare lattice techniques to the perturbative results in this 
regime, because the lattice only determines correlation functions at Euclidean time separa- 
tions; an analytic continuation by the Maximal Entropy Method [10] is required to establish 
dilepton production rates. This method does not work well where the relevant correlation 
function varies rapidly with frequency, which occurs near the light cone. 

The problem at low frequencies is also confused. In one of the pioneering papers of 
the HTL resummation technique, Braaten, Pisarski, and Yuan (BPY) computed the small 
energy, zero momentum dilepton production rate, using HTL resummation to obtain the 
full leading behavior [11]. However, lattice results [8, 12] do not seem to agree with their 
results. Arguably, neither the perturbative result, nor the results from the Bielefeld lattice 
group [8], can be correct, at least in the very low energy limit. The low energy limit of the 
dilepton rate is simply related to the electrical conductivity of the Quark-Gluon Plasma, as 
we will review below. On physical grounds, the electrical conductivity should be nonzero 
and finite. However, the lattice results from the Bielefeld group seem to extrapolate to a 
conductivity which is zero-while the results of BPY extrapolate to a conductivity which is 
infinite. 

This paper will re-analyze the low energy region, using weak coupling techniques. BPY 
found a dilepton production rate dW/d A q ~ q^ 4 , corresponding to a current-current corre- 
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lator spectral weight p(Q) ~ 1/?°, for q° smaller than but formally of order gT. We give 
a clear physical picture of how this scaling behavior arises, and then show that it breaks 
down at the scale q° ~ g 4 T, formally outside the domain of validity of the Braaten, Pisarski 
and Yuan calculation. At very small q°, p/q° takes a shape similar to, though distinct from, 
a Lorentzian, so lim 9 o_^. p/q° is finite, consistent with a finite electrical conductivity. Our 
quantitative results in this region can be considered an extension of existing calculations of 
the electrical conductivity [13, 14]. 

More surprisingly, we find that even in the claimed region of validity of Braaten, Pisarski, 
and Yuan's calculation, q° ~ gT, their result for the dilepton production rate is incomplete; 
for q° <C gT, it is low by about a factor of about 4. This is because the power counting 
underlying their calculation is flawed; there are diagrams which were discarded in their 
analysis which in fact contribute at the same order in the coupling as those they have 
included. We provide a complete calculation for the case q° <C gT. 



II. SPECTRAL WEIGHT, CONDUCTIVITY, DILEPTON RATE 

First let us review the relation between the current- current spectral weight, the conduc- 
tivity, and the dilepton production rate. Dilepton pairs are produced electromagnetically. 
The amplitude arises at second order in the electromagnetic interaction, requiring one inser- 
tion of eA^ Jje pt , with J^ ept the leptonic electromagnetic current, and one insertion of A u Jq CD 
with Jq CD the hadronic electromagnetic current, 

^cz> = I]QwV (2-1) 

(We use [—+++] metric convention.) The differential rate to produce a dilepton pair of 4- 
momentum g M from a QCD plasma described by a density matrix p depends on the squared 
matrix element, 

E / 7 Tr P eJ% CD (q)M-<l) e4 pt (-r')A u (r')\l + (p)l~(q-p)h) 

, J prr 1 ^ 

x(l + (p)r(q-p)h\eJ» CD (-q)A,(q) eJ^ t {r)A u {-r) , (2.2) 



where \h) represents an arbitrary hadronic state. The gauge fields contract to give propaga- 
tors, rj^ u S(q — r)/q 2 , and the leptonic part of the calculation can be evaluated perturbatively 
[3]. The sum on the hadronic final state \h) is the identity. Performing the leptonic part of 
the calculation, the differential dilepton production rate per unit 4-volume is 

where W(q)d 4 q represents the production rate per 4-volume of dileptons of 4-momentum 
q. The required correlation function is called the Wightman current-current correlation 
function, 

KM = J d'xe^Tr pJ^MO) = J d 4 xe^ x (J,(x) J V {G)) P . (2.4) 
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The current correlation is proportional, at lowest order in e and to all orders in g, to the 
one-particle irreducible photon self-energy, e 2 U^ u (q) = ITS (5). 

The spectral weight is defined using the commutator of the current operators, 

p^(q) = J d 4 xe^ x ( [j M (0) , Ux)] )^ = -ft< (5) + fl^(-q) . (2.5) 

In a general density matrix there need not be a simple relation between p^ u (q) and Il < (g), 
but for a thermal density matrix, they are related by the KMS condition [15]: 

PM = (f* - l) n< (g) , or ft< fa) = ^J^jPM = n b (q°)p^(q) , (2.6) 

with rib the Bose statistical function. Note, for q° < 0, this relation is, p^ u {q) = 
— [l+nb(\q°\)]Il< u (q) . Therefore, in terms of the spectral weight, the dilepton production 
rate is 

dW(q) = a 2 p(q) 
d A q ~ 67T 3 g 2 (e^o - 1) ' 1 ' ' 

where p(q) = g^p^{q)- 

Note that the spectral weight p^ v {q) is odd in q, p^q) = — Pjxv{— Also, its relation 
to the (well behaved) Wightman correlator, Eq. (2.6), means that p(q) will vanish linearly 
as q° — > 0. This sometimes makes working in terms of p awkward; it is p^u(q)/q which is 
well behaved near q° = 0. In particular, consider the zero momentum limit of the time-time 
spectral weight, poo(<? = 0, q°). It is related to the Wightman correlator IlQ (<f = 0,q°), 
which is the frequency transform of the real-time Wightman correlation function, 

n<(g = 0,t)= J d 3 xfl< (x,t) = J d 3 x(J (x,t)J (0,0)) . (2.8) 

But this is just the charge susceptibility; writing the spatial volume as V, and using trans- 
lation invariance, 

J d 3 xfl<(x,t) = J d 3 x(J (x,t)J (0,0)} = ^J d 3 xd 3 y(J (x,t)J (y,0)} = i(Q(fjQ(0)) , 

(2.9) 

with Q = f Jod 3 x the total charge. Since charge is conserved, Q(t) = Q(0). The suscepti- 
bility is defined as, 

XQ^^iQ 2 ), (2.10) 

so clearly 

U<(q = 0,t) = XQ. (2.H) 

Fourier transforming, 

^o(Q = ^1 ) = XQ2n5(q°); (2.12) 

the correlator vanishes away from q° = but has a delta function there. The spectral weight 
Poo will also vanish away from zero, but the spectral weight is multiplied by l/rib(q ) = 0. 
However, p/q° retains the delta function and ensures that frequency integrals involving p o 
will correctly be proportional to the charge susceptibility. 
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Now, consider electrical conductivity. It is defined in terms of the current response 
induced by a static electrical field, 



(eJt M ) = a(E), (2.13) 

where (E) is a small, externally imposed electrical field. We assume here that we are in the 
local fluid rest frame. 

The value of the conductivity can be related to the current-current correlator discussed 
earlier via a Kubo relation [16]. Essentially, an external electric field, E^t) = EiRe e~ iq0t , 
is added to the Hamiltonian. This involves introducing a JiAi term, with Ei(t) = doA iy as 
a perturbation to the Hamiltonian; 

H(t)=H + H I , Ej = [ d 3 x eAi (x, t) J, (x, t) , A = , (2. 14) 

J ~W 

where the real part is to be taken at the end. Here H is the full QCD Hamiltonian, 
but neglecting electromagnetic interactions and without any external electrical field. The 
current is evaluated in this perturbed system, fixing retarded boundary conditions such that 
the initial conditions at t = — oo are the equilibrium density matrix. The conductivity is 
given by the ratio of current to electric field in the small frequency limit. Since we are 
interested in weak electrical fields, we may work in the interaction picture and expand to 
first order in the perturbation; 

( Mt = 0))pert = ReTrp th (^exp J -iH^dt^j J^exp J -iH^dt^j 

,o 

= Re / dt'Ti Pth (iH I {t')J l -J l H I {t')) 

J —oo 

= Re J d 3 x J dt 1 '6 (-0 Tr p th -ie J^O, 0) , Aj(x, f)Jj(x, t') .(2.15) 

Since Aj is an external field, we can pull it out, leaving the current-current commutator. 
Also substituting A for E and dividing it through, we find [17], 

^ = 5 Re ng(g = 0,g°), (2.16) 
H/j q 

where n R is the retarded correlator, the one including the O(i) time ordering step function. 
Its real part is related to the spectral weight by 1 

p^(g) = 2Ren R (g). (2.17) 



The spectral weight is the discontinuity of the propagator across the real axis, equivalent to the sum of 
the retarded and advanced propagators. When time translation invariance holds, as in equilibrium, the 
retarded and advanced propagators are related as Ga = —G^, so p = 2 Re Gr. The retarded correlator 
is frequently introduced with a factor of i in the definition, in which case it would be twice the imaginary 
part which equals the spectral weight. 
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Rotational invariance of the plasma ensures that Pij(q = 0) = ^S^pkhiq = 0)> an d so 



3" 

e 2 n 



a=-lim- Pii (q = 0,q o ). (2.18) 
6 g °->-o g u 

This is related to the Wightman correlator introduced previously by 

e 2 .. Pu(q = 0,q°) [3e 2 ~ 

* = T Jte, 7° = "6"fe^ (9 = °' g } • (2 ' 19) 

It is expected that in a real, dissipative system such as QCD at finite coupling, the 
electrical conductivity should be nonzero and finite. Therefore the small q° limit of pa 
should be pa oc q°, and the limit of W(q) should be W(q) oc (g°)~ 2 . 

The weak coupling expansion for such correlation functions near zero frequency turns 
out to be quite involved [18-20]; nevertheless, leading-order weak-coupling results for the 
conductivity of QCD exist [13, 14]. 



III. SPECTRAL WEIGHT AND KINETIC THEORY 

Our emphasis in this paper is on the spectral weight at zero momentum and at frequencies 
q° < gT. The scale gT is relevant because the dominant scattering process in the plasma, 
soft gluon exchange, involves energy and momentum exchanges of order gT (since lower 
exchange momenta are screened). Therefore, the physical time scale for such a scattering is 
~ 1/gT, parametrically shorter than the l/q° time scale separating the current operators in 
the spectral weight. Therefore, to the extent that one is willing to perform an expansion in 
q°/gT <C 1, one can take the scattering events to be instantaneous. In this case, the current 
response can be calculated using kinetic theory. 

The way to use kinetic theory to find the spectral weight is to take the derivation of 
the Kubo relation, discussed above, and use it in reverse; one computes the spectral weight 
by evaluating the current in a theory with an externally imposed electrical field. Because 
the weakly coupled theory has well defined quasiparticles, and because they dominate the 
current, energy density, and so forth, a kinetic theory (Boltzmann equation) description 
should capture the leading order behavior of the current in the electric field background. 

We first introduce a distribution function, f(p,x,t), characterizing the phase space den- 
sity of the quasi-particles. 2 This function satisfies the Boltzmann equation, 



d d d 

di + Vp 'dx~ + cxt 'dp~ 



f(p : x,t) = -C[f], (3.1) 



where v p = dp /dp is the velocity of an excitation with momentum p, F cxt is an external 
force acting on the excitation, and C[f] is the collision integral which characterizes scatter- 
ings with other excitations in the plasma. Since we are primarily concerned with weakly 



2 Our convention is that / is the phase space density for a given spin and color, not summed over spins and 
colors. 
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coupled QCD at small or vanishing quark mass, we can take the velocities to be lightlike, 
v p = p, with m 2 /T 2 and g 2 corrections which we will neglect. 

We are interested in the case where everything is statistically spatially homogeneous, but 
F cxt = QeE is determined by an externally imposed electric field. For this case, the lefthand 
side of the Boltzmann equation can be evaluated, if we linearize in the electric field (which is 
all that is necessary to obtain the current- current correlation function). At linearized order, 
we can write 

f(p, x, t) = f {p) + Sf(p, t) , 5f(p, t) = 5f(p)e-^ 0t , (3.2) 

with 5f(p) in general complex. (When we take the real part of the retarded correlator to 
get the spectral weight, only the real part of Sf will contribute.) Since the electric field is 
already linear order in E, the term containing E need only be expanded to zero order in 
the electric field, 

F- -^/.W-geE- if 

= (3.3) 

where the upper sign is for bosons and the lower sign is for fermions. In fact, since the 
electric charge Q = for gluons, only the lower sign is needed. 
The time derivative acts only on 5f, giving 

d t f(p,t) = -iq°5f(p,t). (3.4) 

Since we know that Sf arises at linear order in the electric field, it must be proportional 
to E. Rotational symmetry completely fixes the possible form of this angular dependence, 
Sf(p) oc E • p 5f(p), with 5f(p) an unknown function of p = \p\ only. It is convenient to 
parametrize this unknown function as 

Sf{p) = /o(l-/o)/i(p) = ^^/o(l-/o)x(p) , (3-5) 

with xip) unknown and to be determined. 

The collision term represents the change in the number of particles of momentum p due 
to scatterings and is a sum of two terms; a loss term represents the scattering of a particle 
with momentum p, out of this momentum state; the gain term accounts for scatterings 
which generate a particle of momentum p. For 2^2 scattering processes, the collision 
term is 



C 2 r 2 [f](p) = jr^Y. I \M%(p,k;p',k')\ 2 (2ir)W\P + K-P'-K') (3.6) 

'n P )f{k)[i ± /v)][i ± f d m] - np')f d {k')[i ± ± f\k)\) , 



where the first product of statistical functions is the loss term and the second is the gain 
term. 3 Here |-M^| 2 is the spin and color summed squared matrix element for ab — > cd 



3 We write the collision term as loss minus gain, and put a — sign in front of it in Eq. (3.1), so that its 
linearization will be a positive operator. 
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and J k = J d 3 k / (27r) 3 2k° is the Lorentz invariant phase space integral. The leading factor 
l/4p°z/ a is the usual l/2p°, a \ jv a to replace the initial spin and color summed |.M| 2 with 
the spin and color averaged one, and a 1/2 to eliminate double counting final states, cd — > cfc 
(or as a final state symmetry factor when c — d). 

Just as we have linearized the lefthand side of the Boltzmann equation, now we must 
linearize the collision term. Because of our particularly convenient choice Sf = / (1 ± fo)fi, 
the loss minus gain difference of population functions takes a very simple form, 



r(p)f b (k)[i±f c (p f )][i±f d (k f )} - n P f )f d (k')[i±n P )}[i±f b (k)} 

= /o a b)/o(fc)[i±/o(pO][i±/o d (fcO]{/i(p)+/i(fc)-A(p / )-/i(fc / )} • 

Substituting 

/i(fc) = Qe^x(k) 

looks like it will make the result complicated; on the left side of the Boltzmann equation, only 
E ■ p appears. Fortunately, the phase space, matrix element, and equilibrium distribution 
functions are all invariant on rotations about the p axis. Averaging over such rotations 
replaces k with (k ■ p)p, so after such averaging, the Boltzmann equation becomes 



(3.7) 



(3.8) 



QT/ (1 - f )p + iq°Qf (l - f )px(p) = 



P 



E 



Mf d ( P ,k;p f ,k 



bed Jk P' k ' 

x(2n)^\P + K-P'- K')f a (p)f (k)[l ± / c (p')][l ± f$(k')] 
xp • \QPX(P) + Q b k X (k) - QcP'xip') - Qdk' X (k') 



(3.9) 



or 



QT = -tq°Q X (p) + 



1 



-E 



bed " k P' k ' 



/o(l-/o)4p^ u 

x/o 8 (p)/o 6 (fc)[l±/oV)][l±/o(fe')] 

X 



M? d (p,k;p', k')\ 2 (27r) 4 5 {4 \P+K-P'-K') 

(3.10) 



Qxip) + p ■ kQ bX (k) - p ■ p'QeXip') ~ P ■ k'QdX(k') 
Symbolically, this can be expressed as, 

S(p)=[-tq° + c]x(p), (3.11) 

where C, the linearized collision operator defined above, is a nonlocal operator in p space. 



This can be formally solved for x by inversion, x(p) 



-iq 



C]-'S(p). 



We will solve this (integral) equation in subsequent sections. First, though, we show what 
the solution has to do with the current-current correlation function which we are after. The 
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current, in the kinetic theory approximation, is 



Ji = 




(3.12) 



In passing to the final line we used that the angular average of pipj is 5^/3. We see 
as expected that j oc E; (the integral of) the function x expresses the coefficient of this 
proportionality. Since x is i n general complex (see Eq. (3.11)), the current typically lags the 
electric field; the dot product j ■ E/E 2 depends on the real part of x{p)- Using Eq. (2.16) 
and Eq. (2.17) from the last section, we see that the spectral weight is 



Therefore the finite frequency solution to the Boltzmann equation yields the spectral weight 
at finite frequency. This is a generalization of the technique used to determine the electrical 
conductivity in [13]. 

IV. DIAGRAMS, PHOTON EMISSION, AND BOLTZMANN EQUATION 

Before computing the spectral weight, we pause to discuss the relation between what we 
have just done and the usual diagrammatic approach to the photon production rate. In 
doing this, we will see the physical origin of the flaw in the treatment of BPY. 

We begin by emphasizing the range of validity of each calculation. The Boltzmann 
calculation is most naturally phrased in terms of frequencies q° ~ g A T, since this is where 
the collision and iq° terms in the Boltzmann equation, Eq. (3.11), are of the same order. 
However, the treatment only breaks down (at weak coupling) when the "instantaneous 
scattering" approximation in the Boltzmann treatment ceases to be valid, so it is valid 
provided q° <C gT '. 

The treatment of BPY was billed as a complete calculation for q° ~ gT. As such, it 
should be valid over a parametrically wide range. By treating q° as formally O(gT), the 
authors' approximations could become uncontrolled when q° differs from this estimate by a 
power of g, for instance, for q° ~ T or q° ~ g 2 T. However, the region g 2 T <C q° <C gT should 
be within the validity of their approximations. The two calculations overlap in this range, 
and the absence of any change in behavior at the g 2 T scale in the Boltzmann treatment 
shows that the BPY calculation should in fact be valid for any scale q° 3> g 4 T. However, 
their treatment certainly fails below this scale (and they never claimed otherwise). Hence, 
their treatment is valid for g 4 T <C q° <C T. Therefore, both treatments should be valid in 
the range g 4 T <C q° <C gT, and we can compare with their results in this regime. 



Pij(q°) = 25y 
o° ~ 3T 



5>a<£ / / (l-/o)Xa(p) , 



S = 2 5>$ / /o(Wo)Xa(p) ■ (3.13) 
3 _ J J v 
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FIG. 1: The amplitude for a Compton scattering to produce a photon involves a sum over where 
the external photon attaches. 



Because q° ^> C in this regime, we can simplify the Boltzmann treatment by making an 
expansion in large q° in Eq. (3.11): 



-iq° + C 



i 1 A 
q° {q ) 2 



(4.1) 



In other words, we can treat iq°x as the dominant contribution in Eq. (3.11) and solve the 
equation iteratively. The leading term, \ = iT/q°, is pure imaginary and does not contribute 
to the spectral weight. The next contribution arises by substituting the leading contribution 
into the collision term, and gives 



QaXa 



(2) 



■ / £ \M 2 \{2^5\P+K-P'-K')f{ P )f{k) 

Jp'kk' 



(Wo(Wo)4piW«« M 

x [l±f(p')] [l±f(k')] Q a + p-kQ b -p- p'Qc - P • k'Qd 



(4.2) 



which, substituted into Eq. (3.13), yields a spectral weight of 



P 



(q°) 2 T 



(q°) 2 T 



PP'kk' 2 ahri . 



~ E \Mt\\^f5\P+K-P'-K')f(p)f(k)(l±f(p')(l±f(k')) 



XQaP ■ [QaP + Qbk - Q c p - Q d k' 



pp'kk' 4 abcd 



7 E l-M c 1| 2 (2vr) 4 5 4 (P+K-P'-A-)/(p)/(fc)(l±/(p / )(l±/(^)) 



x (QaP + Qbk - Q c p' - Q d k' 



(4.3) 



where the last expression is a symmetrization of the previous one over the indices abed. 

This expression has a simple diagrammatic interpretation. The momentum integration, 
particle sum, matrix element, momentum conserving delta function, and population function 
give the total rate per 4- volume for a scattering process to occur. The factor of 1/4 corrects 
for overcounting a <-> b and c <-> d (or is a symmetry factor if these labels are the same). 
A factor of e 2 times the l/(g ) 2 times the final factor represent the probability that the 
scattering process should lead to an off-shell photon emission with frequency q°. 

To see this, consider how a typical scattering rate is modified by the addition of a soft 
photon. For concreteness consider Compton scattering. The photon can attach to any 
charged external leg or internal propagator, see Figure 1. However, it is only the amplitude 
when it attaches to an external leg which is l/q° enhanced, as we now discuss. 

Consider the diagram where a soft photon is emitted from one of the quark external legs. 
For this discussion we write the photon momentum as R to avoid confusion with the charge. 
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Take the external state momentum to be P, so the momentum on the internal propagator 
is L = P =1= R (where =F is — for initial and + for final state particles.) Because the particle 
is on shell, P 2 = 0. 

Writing the polarization vector for the photon as e, the diagram with the photon emission 
on an initial state leg differs from the diagram without any emission by the substitution 

u(P) - ^ [-Qi\ u{P) (4.4) 

and in a shift to the kinematics of the rest of the process, P — > L. Because R is small, the 
change to the kinematics of the rest of the diagram represent a small correction and can be 
neglected. 

Since L = P - R, 



L 2 P 2 + R 2 - 2P ■ R ~ -IP ■ R 2p°r° 



(4.5) 



Anticommuting f with ft and using that fu(P) = (the Dirac equation), Eq. (4.4) is 
approximately 

<P)^%^<P) = ^QP'*<P)- (4-6) 

The same thing happens when we consider an outgoing external leg, except that it is +2P-R 
rather than —2P ■ R which appears, so the overall sign is opposite. The term where the 
photon attaches to an internal leg is proportional to the virtuality of that leg, which is at 
least O(gT) and is therefore subdominant at small frequency. Therefore the amplitude for 
the bremsstrahlung process is identical to the amplitude for the regular 2 <-> 2 diagram, but 
with an overall multiplicative factor of Sin-out P ' e - On squaring an d summing over 
polarization states (and there are three, since the photon is off-shell and in the pure time 
direction), we recover the factor e 2 /(q ) 2 x (Sin-out Qp) 2 which we found above. 

Therefore the result, Eq. (4.3), is just what we would obtain diagrammatically by con- 
sidering bremsstrahlung processes, after accounting for the fact that the photon kinematics 
are soft. 

With the exception of Coulomb scattering, all 2 <-> 2 processes in the plasma have a rate 
~ a 2 T A up to logs. Coulomb scattering is faster by a factor of l/a s , but only because of 
small angle processes for which (p — p') 2 ~ a s ; therefore its contribution to Eq. (4.3) is 
only logarithmically larger than a typical 2^2 process. Namely, each such process gives a 
contribution of order 

p a 2 T 3 dW a 2 a 2 T 3 



(4.7) 



q° (q ) 2 ' d 4 q (g ) 4 ' 

which is parametrically 0(a 2 /T) for q° ~ gT. 

Our result is the same order, and shows the same q° dependence, as the term found 
by BPY. However, the diagrams they include in their analysis only give rise to photon 
production from Compton scattering and its crossings, not, for instance, from Coulomb 
scattering. Therefore their calculation must be incomplete. 

Let us make a more detailed comparison. The complete result of BPY, for general quark 
number and charge, is [11] 
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dW(q) 
d 4 q 



^0 J k2dk J duJ J du'n f {uj)n f {uj')5{q - uj - u') 
x<! s|l n 1 p+(u,k)p-(u',k) 



p=0 



where 



2kq 



u 2 + J 2 - 2k 2 - 2m 2 
2kq® 

uj 2 + u' 2 - 2k 2 - 2m 2 \ 2 



2kq° 



p + (uj,k)p+(uj',k) 
p-(uj,k)p-(uj',k) 



x 



p±(uj,k) 



u 2 -k 2 
2m\ 



[6(u - u±(h)) + 5{uo + uo T {k))} + /3±(u, k)Q(k 2 - uo 2 ) , 



(4.9) 



and 



P±(u,k) 



(m 2 q /2k)(lTuj/k) 

{uo T k - (m 2 /k)[Q (u/k) T Qi(u/k)}} 2 + [(nm 2 g /2k)(l =F u/k)]< 



(4.10) 



We are interested in the small q° behavior, for which uj 1 ~ — uj. In this case the prefactors on 
the spectral weights contain large l/(q ) 2 enhancements, but only when uj 2 ^ k 2 , meaning 
that only the f3±, or cut, piece of p± contributes appreciably. In this case the integrand 
behaves as m^/A; 4 ^, indicating a large k log divergence, which BPY remarkably did not 
comment on. This is the familiar log divergence due to the s/t nature of the Compton 
cross-section, seen from the infrared side. 

Attempting only to compute the coefficient of this log divergence, we may expand in 
\uj\,k 3> m q . (Here m 2 = [C-p = A/3]g 2 T 2 /8 is the effective thermal mass of a quark.) The 
expansion is relatively straightforward and leads to 



dW 
~¥q~ ( 



0) 



2a 2 ZQ 2 3 ™ 



n 3 q 2 



16nq 2 





r k dcu 


fi--l 


1 dk 


l_ k T 


k 2 _ 



(4.11) 



The corresponding log divergent coefficient in the spectral weight is 



p(q = 0,q°) 



x 



4^T 



(4.12) 



where we have taken the upper and lower limits of the k integration to be T (where popula- 
tion functions enter) and m q (where the full expression is required and the result is smaller). 
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It is not clear to us how this divergence was handled in the numerical work presented by 
BPY [11]. 

Now we approach the same calculation using Eq. (4.3). We will show as claimed that 
the BPY result coincides with the contribution from Compton type processes. The relevant 
summed matrix element is 12J2Q 2 x 8Cfs/t, with s, t the usual Mandelstam variables. 
Here 12 is the number of degrees of freedom per quark flavor (2 spins times 3 colors times 
particle/antiparticle) and 8Cps/t is the matrix element squared summed on all other species' 
spins and summed over Compton and annihilation processes. The phase space integrations 
can be reduced using the technique of [14], to give 

P _ 12 E Q 2 x8Ci f dk f k du rd A f dp f dp2 S 



g° 2%s g 2 T J Q J k J Q 2 Wo Jo t 

X//(Pi)[1+/6(Pi)]//(P2)[1+/ 6 (P2)] X (p 1 -p 2 f, (4.13) 

where we have already made the approximation u>, k <C T in writing the limits of p integra- 
tions, population functions, and the vectorial term. Physically, k, cu represent the exchange 
momentum and energy in a Compton type process and Pi,P2 represent the energies of the 
two participating (anti) quarks. The Mandelstam variable ratio and the required angle are 
given approximately by 

- ~ -p-(l-cos0), p 1 .p 2 = l-__~l __(i_ cos 0). (4.14) 

The pi, p 2 , and integrations are then completely straightforward, and they yield 



1 '' •"• i' 1 ,, 

j.2 ' ( 4 - 15 ) 



P_ _ io\"n2 3m t fdk [ k du [ _ u 2 

f~ l2 ^ Q itotfrj k J_ k k - 



which is identical to Eq. (4.11) and leads to the same result for the spectral weight as 
Eq. (4.12). 

We see that the leading small q° behavior in the BPY calculation corresponds to evalu- 
ating just the Compton-type processes in the kinetic theory approach. Their calculation is 
therefore obviously incomplete. What went wrong in the BPY calculation? Their mistake 
was to assume that the loopwise expansion would be valid for photon production, modulo 
HTL corrections. This would be true were they considering a hard (q° ~ T) virtual external 
photon-in this case the first corrections are known and are 0(a s ) [7]. However, considering 
the q° ~ gT regime, low-loop calculations can be suppressed by kinematic restrictions which 
disappear at higher loop orders. For instance, consider the 1-loop diagram analyzed by BPY, 
a two-loop diagram, and an example of the 3-loop diagrams which we find to dominate at 
small momentum, all shown in Fig. 2. Besides the explicit e 2 in the one-loop diagram, the 
fermion lines are kinematically constrained to carry momentum ~ q° ~ gT so that they 
may both be simultaneously on-shell; since the photon line attaches to this low-momentum 
line, the diagram has an extra factor of g 2 ~ a s and is 0(aa s ). This is not enhanced by any 
Bose stimulation functions because the soft lines are fermionic. 

For the three- loop diagram, on the other hand, there is an explicit factor of e 2 g 4 , but 
the softness of the incoming momentum makes no restrictions on the sizes of the three 
loop momenta, which can each be 0(T) and at large angles with respect to each other 
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FIG. 2: One loop diagram considered by BPY, a two-loop diagram, and an example three-loop 
diagram which we find contributes at the same order for q° ~ gT. The hashed diagonal line 
represents a cut; lines crossing the cut are put on-shell. 

without preventing the four cut lines from being simultaneously on-shell. The momentum 
factors therefore lead to a T 2 behavior. Furthermore, there are two uncut, almost on-shell 
propagators which lead to a 1/(P • Q) 2 ~ l/g 2 enhancement, giving a rate which is 0(aa s ), 
the same order as the 1-loop diagram. This is true even if all of the lines in the 3-loop 
diagram are hard-the need to add loops is not just the well-understood need to resum 
HTL's on soft lines. 

Almost the same is true of the 2-loop diagram, except that if the lines are all hard then 
there is a collinear restriction for the three cut lines to be simultaneously on-shell, reducing 
the phase space by 0(a s ) and again making the rate 0(aa s ). Higher loop diagrams are 
suppressed with respect to these, however, because loop factors continue to appear but no 
additional kinematic constraints are lifted, since there are no remaining kinematic constraints 
on the 3-loop diagrams. This argument fails when q° ~ g 4 T, at which point a large class of 
diagrams begin to contribute (which we effectively resummed in using kinetic theory). 

V. KINETIC THEORY EVALUATION OF SPECTRAL WEIGHT 

We are now ready to proceed with the evaluation of Eq. (3.11) and Eq. (3.13). This 
involves solving a non-trivial integral equation, which cannot be handled analytically. We 
must revert to a numerical treatment, which can however be made highly accurate by using 
the variational techniques applied to conductivity in [13, 14]. 

The quantities of interest in the evaluation of Eq. (3.9), T/ (l— fo)p = S(p) and pxip), 
are general vector functions of p (with particularly simple angular dependence). The colli- 
sion integral and q°fo(l— fo) can be considered as linear operators on such functions (with 
q°fo(l— fo) being an almost trivial operator). It is useful to consider such functions of p as 
elements in the vector space of all functions of p with the same angular dependence, and to 
define an inner product on such functions / and g, namely 



The utility of this inner product is that the collision operator is real, positive, and symmetric 
under this inner product, and the quantity we want, Eq. (3.13), is the inner product of two 





14 



functions, (^x,S^j. Therefore, a variational approach for determining x(p) wm have errors 
which are quadratic in the quality of the variational Ansatz. 
With this in mind, we define the functional 



Q[x] = (X,S) -\{x, [C-tq°f (l-f )] x) 



(5.2) 



Extremizing Q[x] with respect to x precisely yields the linearized Boltzmann equation (3.11), 
and the extremal value of Q is simply related to the quantity we want; 

2 



P_ 

q° 



Re ( x , S^j = 4 Re Q 



cxtr • 



(5.3) 



Unfortunately, the matrix C — iq° fo{l— fo) is not positive definite, so the real part of the 
extremal value is not necessarily a lower bound for the result at finite q°. 

The extremization of the functional Q[x] may be performed in several ways. We will 
follow [13, 14] and treat the problem variationally. For an exact solution, we would need to 
work with an infinite dimensional space of arbitrary functions xip)- However, by choosing 
a suitable set of basis functions, (j)^ m \p), and performing the maximization in this finite 
dimensional variational vector space, we are able to get a highly accurate approximate 
result. We will consider the Ansatz, 



K 



(5.4) 



m=l 



where K is the size of the basis set considered and Xm are the variational parameters used 
to tune the maximization of Q[x\- Inserting our Ansatz into the source and collision terms 
of the functional yields linear and quadratic combinations of the arbitrary coefficients Xw 



x,s 
xXx 
(x,?7o(i-/o)x) 



EY a S a 

a,m 

y a c ab y b 

/ j / j Am mnAn > 
a, ra b,n 

/ j / j A.mqmn u abXr. 



(5.5) 
(5.6) 
(5.7) 



In these expressions, 



S a , (m) 



q° 



and 



a,m b,n 

g°^ m ),/o(l-/o)0 {n) / 

Restricted to this subspace, and treating x, S as vectors and C and q° 
as matrices with index (a,m), the functional becomes 

Q[x]=x v ~s-W\c 



if\x 



(5.8) 



We note that the previous equation has both a real and imaginary part to it. Since we 
ultimately want the real part, it will be convenient to write the complex structure in matrix 
form; 



Q[x\ 



S,0 



Re x 
Im x 



Re x , Im x 



C q° 



Re x 
Im x 



(5.9) 
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Finding the approximate form of the spectral weight (the real part of the extremum) is 
now an elementary linear algebra exercise and is given by 



Q max = lS T M- 1 S, 



where, 



S = 







M = 



UJ r . 



(5.10) 



(5.11) 



The next step is to choose suitable trial functions. As discussed in [14], a good choice for 
the basis functions will consists of 



(m) 



(P) = 



ip/T) 



m—l 



'l+p/T) K - 2 ' 



m = 1, K. 



(5.12) 



This set spans both the function S itself, and the variational Ansatz, x{p) V-, which 
has been found to work remarkably well for conductivity [21]. These functions are not 
orthogonal, but they do not need to be. At any finite K they are not complete, but they 
span nested subspaces (the space spanned at one K sits in the space spanned at a larger 
K), and as K — ■> oo, the basis in fact becomes complete. The integrations to establish 
S n can be performed as straightforward quadratures integrations; and while the integrals 
for C mn are less straightforward, they can be reduced to 4-dimensional integrals and also 
carried out by quadratures. One complication is that even after scaling out the explicit 
factor of g 4 from the collision term C, it retains a logarithmic dependence on g because 
Compton and Coulombic scattering processes give infrared logarithmic divergences, cut off 
by plasma screening; and the plasma screening cutoff scale varies with g. Both our approach 
to the numerical integrations, and the details of how screening cut off IR singular scattering 
processes are dealt with at length in [14], so we will not present the details again here; the 
main new feature is the complex matrix structure which we have already clarified. 



VI. RESULTS 

The scale where the spectral weight shows nontrivial behavior is q° ~ g 4 T, since this is 
where the collision term, which is explicitly 0(g 4 ), becomes comparable to q° in evaluating 
Eq. (5.10). In this case p/q° ~ T/g 4 . The weak-coupling behavior of the spectral weight, 
scaling out this explicit behavior, is displayed in Fig. 3. The result depends on the number 
and charge of quark flavors, and it also retains some dependence on the coupling, because of 
the screening masses mentioned above. The simplest functional form with the same behavior 
as p/g°-an even function with a finite value at zero and I/c/q tails-is the Lorentzian. Though 
we have not shown it, the large mo (large coupling) curves in Fig. 3 are quite well fitted by 
a Lorentzian, but the small curves are less so. 

We also present the functional form of p/q° as a function of q°/T for 3-flavor QCD in Fig. 
4. The figure shows the small q° behavior we have found for three values of the coupling 
a s = 0.1,0.2,0.3, together with the leading-order large q° behavior. The region where the 
curves cross is parametrically gT; neither calculation works here. The calculation of BPY 
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Photon energy q°/afl 



FIG. 3: Spectral weight vs. q with the dominant a s dependence factored out. For Nf = 2, the 
two values of mj) correspond to a s = 0.06 and 0.24; for Nf = 3 they correspond to a s = 0.05 and 
0.21. 



o 10 



U 

-t-J 
a 

0) 

a 

CO 



5 - 




left: spectral peak at 
a =0.1, 0.2, 0.3 
(sharpest to broadest) 

Below: leading large q° 
(Born term) 



2 3 4 5 

Energy q°/T 

FIG. 4: Dominant large q° and small q° contributions to the spectral weight. The rising curve is 
the leading-order Born term; the three curves sharply peaked at the origin are our results for a 
3-flavor plasma with a s = 0.1,0.2,0.3. As the coupling is reduced, the peak at small q° becomes 
taller and narrower, but the area underneath is unchanged. The Born result is valid where it is 
much larger than the peak, and the peak is valid where much larger than the Born term. In the 
overlap region (where q° ~ gT) the detailed behavior is not known. 
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FIG. 5: Large q° asymptotic of the kinetic treatment (which is the small q° asymptotic of the 
q° ~ gT treatment) as a function of the gauge coupling, for Nf = 2, 3. 



was intended to establish p/q° here, but we have already shown that it is incomplete, so at 
present there are no leading-order perturbative results in this region. 

As we already emphasized, it is especially easy to extract the q° 3> g 4 T asymptotic of 
the above curves, which is supposed to agree with the q° <C gT asymptotic of BPY. We 
obtained an explicit expression in Eq. (4.3). We present this asymptotic as a function of a s 
in Fig. 5. It can also be expanded in the logarithm of the coupling; since it is the collision 
operator and not its inverse which appears, this expansion terminates rather than containing 
all orders in inverse logs, as is the case for the conductivity [14]. Explicitly, 



4 [ 3 'T « 9 « « 9 T] = 



4tt 16(6 + JV f 



C 



3 3tt 

-0.010 N f = 2 
-0.039 Nf = 3 ' 



ln- + C) , (6.1 ) 

9 



The first term in the square bracket arises from Compton-type processes, the second is 
from Coulombic scattering; they were obtained by evaluating Eq. (4.3) using the techniques 
of section IV of [13]. We saw in Eq. (4.12) that the treatment of BPY reproduces the 
first but not the second term. Unfortunately, the range of validity of this expansion is 
very narrow, since it becomes negative before a s = 0.1; the leading-log treatment does 
particularly poorly here. Physically this is because the dominant momentum of a charge 
carrier is small, and at physically relevant couplings the typical scattering process changes 
its direction substantially. The case is worse than for transport phenomena because they 
are dominated by larger momentum particles. 

Fig. 4 shows that the spectral weight has a peak near q° = 0. As the coupling increases, 
the peak spreads out but becomes lower. But even at quite large coupling, a s = 0.3, the 
peak is still very close to zero frequency. Further, the area under the peak remains the same 
regardless of the value of a 8 , because the small q° result we have obtained obeys a sum rule. 
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To see this, consider the integral over all q° of p/q°; 

Jdq°-^ = J dq° Re 2(5, [C - iq° / (1- fa)}' 1 s) . (6.2) 

Inserting a complete set of eigenfunctions A of the operator C, with eigenvalue A, this integral 
becomes 

/V 4 = £ / ^° 2Re ( 5 > P-wVoCi-zo)]-^) (a, 5) 

9 A 

= £ / dq° 2 Re (5 , [A - ^(l-/,,)]-^) (a , s) . (6.3) 

A ^ 

The operator C has been replaced by a number A. Therefore, inside the integral implicit in the 
inner product, the inverse, [A — iq° fo{^~ fo)] 1 , is an ordinary (rather than operator) inverse, 
and there is no obstacle to performing the q° integration. The integral is straightforward 
and yields 

/ *° £ = E (* • A) (A .S) = 2,(S , . (6,, 

Explicitly, the remaining integral is 

VII. DISCUSSION 

We have presented a leading-order calculation of the electromagnetic current-current 
spectral weight, and therefore of the dilepton production rate, in weakly coupled, thermal 
QCD and for dilepton energies q° <C gT. In this regime, the computation of the spectral 
weight requires applying kinetic theory. There is a simple interpretation of the kinetic theory 
result for q° ^> g 4 T in terms of bremsstrahlung from scattering processes in the plasma. 

Our results show that the asymptotic behavior found by Braaten, Pisarski, and Yuan, 
p/q° oc q^ 2 for q° <C gT, is correct down to the scale q° ~ g A T, though the coefficient 
which they found for this behavior is in error. This error was caused by an insufficiently 
careful power-counting; a wider class of diagrams than they considered contribute at leading 
order, because the one-loop diagram has kinematic constraints which disappear at higher 
loop order. 

At the scale q° ~ g 4 T, the spectral weight's behavior moderates, so that lim^o 
is finite, as is required by the finiteness of the electrical conductivity. The weakly coupled 
spectral weight shows a large peak at small q°, but the area under this peak is finite. 

Furthermore, there is a sum rule satisfied by p/q°; its q° integral over the aforementioned 
peak is independent both of the coupling and of the details of the scattering processes which 
determine its detailed form. The significance of this sum rule is that it makes it very hard to 
use Euclidean correlation functions to reconstruct the peak value, liny ^ p(<?°) / Q° — 6o"/e 2 
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(with a the electrical conductivity). Since this is precisely what must be done to find 
the conductivity from lattice gauge theory calculations, we believe that (at least at weak 
coupling) such a reconstruction is extremely challenging. 

To clarify this last point, we refer to the discussion of Aarts and Resco [22]. The Euclidean 
correlator at zero spatial momentum and finite Euclidean time r is determined in terms of 
the spectral weight by an integral, 



The point is that K(t, q°) is essentially flat and r independent near q° = 0, which means 
that G E at any r essentially captures the area under the small q° peak. But this area is fixed 
by the sum rule and carries no information about the detailed shape of the peak. Hopefully, 
at very large coupling the peak becomes broad enough (or even disappears [23]) that the 
weak q° dependence in K(t, q°) is enough to see its structure-though even at a s = 0.3, the 
edge of the range where we believe perturbative methods, the peak remains very narrow, 
with most of its support at q° < 0.5T. So long as the peak is narrow, fitting will be difficult, 
especially without an Ansatz for the shape of the peak. Fortunately, as we mentioned, the 
peaks shown in Fig. 3 for large coupling are surprisingly well fit by a Lorentzian, exactly 
the fitting form advocated by Aarts and Resco [22]. 
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